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ABSTRACT 



Context. Knowing accurate lead abundances of metal-poor stars provides constraints on the Pb production mechanisms in the early 
Galaxy. Accurately deriving thorium abundances permits a nucleo-chronometric age determination of the star, 
j Aims. We aim to improve the calculation of the Pb i and Th ii lines in stellar atmospheres based on non-local thermodynamic equilib- 

, rium (non-UTE) line formation, and to evaluate the influence of departures from LTE on Pb and Th abundance determinations through 

\^ • a range of stellar parameters with variations of the metallicity from the solar value down to [Fe/H] = -3. 

Methods. Comprehensive model atoms for Pb i and Th n are presented. We describe calculations of the Pb i energy levels and oscilla- 
tor strengths. 

Results. The main non-LTE mechanism for Pbi is the ultraviolet overionization. Non-LTE leads to systematically depleted total 
^ ■ absorption in the Pbi lines and accordingly, positive abundance corrections. The departures from LTE grow with decreasing metal- 

■ licity. Non-LTE removes the discrepancy between the solar photosphere and the meteoritic lead abundance. Using the semi-empirical 
I Holweger & Miiller (1974) model atmosphere, we determined the lead non-LTE abundance for the Sun to be logEpbg = 2.09. We 
, revised the Pb and Eu abundances of the two strongly r-process enhanced stars CS 31082-001 and HE 1523-0901 and the Roederer et 

al. (2010) metal-poor stellar sample. Our new results provide strong evidence for universal Pb-to-Eu relative r-process yields during 
course of the evolution of the Galaxy. The stars in the -2.3 < [Fe/H] < -1.4 metallicity range have, on average, 0.51 dex higher 
Pb/Eu abundance ratios compared with that of the strongly r-process enhanced stars. We conclude that the s-process production of 
^ , lead started as early as the time when Galactic metallicity had grown to [Fe/H] = -2.3. The average Pb/Eu abundance ratio of the 

■ mildly metal-poor stars, with -1.4 < [Fe/H] < -0.59, is very close to the corresponding Solar System value, in line with the theoreti- 
cs ■ cal predictions of Travaglio et al. (2001) that AGB stars with [Fe/H] ^ - 1 provided the largest contribution to the solar abundance of 

s-nuclei of lead. The departures from LTE for Thii are caused by the pumping transitions from the low-excitation levels, with £exc < 
<n: 1 eV. Non-LTE leads to weakened Thii lines and positive abundance corrections. Overall, the abundance correction does not exceed 

. 0.2 dex when collisions with H i atoms are taken into account in statistical equilibrium calculations. 

. Key words. Line: formation - Nuclear reactions, nucleosynthesis, abundances - Sun: abundances - Stars: abundances - Stars: 

atmospheres 

ITJi [ 1. Introduction The knowledge of stellar Pb abundances provides a bet- 

• ter understanding of the production mechanisms of lead during 

■ Lead (Z = 82), thorium (Z = 90), and uranium (Z = 92) are the evolution of the Galaxy. Lead is produced by two neutron- 
- ■ ■ the three heaviest elements observed in metal-poor (MP) stars, capture processes on different timescales, i.e., in the slow (s) pro- 

with [Fe/HU < -1. Only one de tection of bismuth (Z = 83) has cess occuring during thermally pulsing asymptotic giant branch 
ever been made in an MP star divans et al.ll2005h. because the (AGB) phase of intermed iate-mass (2-8 Mg) stars (see for ex- 
only suitable line, Bii 3067 A, lies in a very crowded spectral ample. iGallino et al.|[T998l) and the rapid (r) process. The astro- 
region in the near-UV that is difficult to access. All nuclei be- physical site at which the r-process operates has not yet been 
tween ^'^Bi and -'^''pu and which exist in nature, are radioactive, identified. Given that the heavy elements, such as Ba and Eu, are 
Among them, only Th and U have long-lived isotopes with half- detected in very metal-poor (old) stars, supernovae with progen- 
lives longer than 1 Gyr, ti/2C^^Th) ^ 14 Gyr and ti/iC^^V) - itors of 8-10Mq are the most promising candidates forr-process 
4.5 Gyr, and can be detected in spectra of MP stars. enrichment i n the early Galaxy. See, for example, the pioneer- 

ing rev iew of lHillebran dt ( 19781) and also lCowan & ThielemannI 

(120041) for further discussion on the r-process. Theoretical stud- 

Send offprint requests to: L. Mashonkina; e-mail: lima@inasan.ru jj^g ^.process showed that in low-metallicity environment. 

In the classical notation, where [X/H] = log(yVx/WH)«» - heavy s-nuclei like those of Pb are more efficiently produced 
log(A'x/A'H)s™. ^ ^ ^ 
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compared to lighter nuclei due to the higher ratio of free neu- 
trons to Fe-peak seeds, and also due to an extended period of 
s-pro cess nucleosynthesis (Gallino et al. 1998, Travaglio et alJ 
l200lh . It was predicted that the AGB stars with [Fe/H] ^ -1 
made the greatest con tribution to the solar ab undance of s-nuclei 
of lead. According to lTravaglio et alj (1200 Ih . the s-process frac- 
tion of the solar Pb amounts to 81 %. 

The larg est data set on lead a bundances of MP stars was 
compiled by iRoederer et al.l ( l2010h . They found that the stellar 
Pb/Eu abundance ratios form a plateau in the metallicity range 
-2.3 < [Fe/H] < -1.4, and Pb/Eu grows at higher metallic- 
ity. The increase of Pb/Eu at [Fe/H] > - 1 .4 suggests a grow- 
ing contribution to the stellar lead abundance from the s-process 
opening in AGB stars. iRoederer et al.l (1201 Ol) concluded that s- 
process material was not widely dispersed until the Galactic 
metallicity grew considerably, as high as [Fe/H] - -1.4. This 
also means that lead in stars with lower metallicity should be 
of pure r-process origin. It is worth noting t hat despite the 
large s catter in the Pb/Eu abundance ratios in the IRoederer et alJ 
(1201 Ol) stars with -2.3 < [Fe/H] < -1.4, the ratios are, on aver- 
age, 0.85 dex higher compared with the correspondi ng value in 
the st rongly r-process enh anced star CS 3 1082-001 dPlez et alJ 
l2004h . With [Eu/Fe] = 1.6 dHill et al.l2002h . this star experienced 
enrichment from, probably, a single r-pro cess source. Following 
the suggestion of IChristlieb et al 1 (12004 . we hereafter refer to 
stars having [Eu/Fe] > +1 and [Ba/Eu] < as r-II stars. The 
heavy element abundance ratios of each r-II star are expected to 
reflect the relative yields of the r-process. What can be the rea- 
son of the large discrepancy in Pb/Eu between CS 31082-001 
and the -2.3 < [Fe/H] < -1.4 stars? 

In the visual spectra of MP stars, lead is represented by the 
lines of Pb i. Due to a relatively low ionization energy of 7.43 eV, 
lead is mostly ionized in the stellar atmospheres with effective 
temperature hotter than Tgff - 4000 K, and the number density of 
neutral lead can easily deviate from the thermodynamic equilib- 
rium (TE) population, owing to deviations of the mean intensity 
of ionizing radiation from the Planck function. All previous stel- 
lar Pb analyses have been made under the local thermodynamic 
equilibrium (LTE) assumption, and the obtained results should 
be checked for any departures from LTE. 

The detection of thorium permits, in principle, a nucleo- 
chronometric age estimate of the star by comparing the observed 
Th-to-stable neutron-capture element-abundance ratios with the 
corresponding initial values provided by the r-process event just 
prior to the time when the star was born (see, for ex ample, 
LSneden et al.llT996llCayrel e"tani200lLlFtebel et alJl2007h . Given 
that the half-Ufe of ^^^Th is long, stellar thorium has to be mea- 
sured very precisely, with abundance errors of no more than 
0.05 dex, to provide less than 2.3 Gyr uncertainty in the star's 
age. 

This study aimed to improve the calculation of the Pb i and 
Th II lines in stellar atmospheres based on the non-local thermo- 
dynamic equilibrium (non-LTE) line formation and to evaluate 
systematic abundance errors caused by the simplified LTE line 
formation treatment. We constructed, for the first time, model 
atoms for Pb i and Th ii. The non-LTE method and the departures 
from LTE for Pb i are described in Sect.|2l Our theoretical results 
were used to determine the solar l ead abundance (Sect.[3l l and 
to revise the Pb abundances of the IRoederer et alJ (1201 Oh stel- 
lar sample (Sect.gl l. Having added the r-II star HE 1523-0901 
(iFrebel et al.ll2007 ). we found that the Pb/Eu abundance ratios 
of the two r-II stars are well re produced by the waitin g-point r- 
process model as presented by IRoederer et alj (l2009l) . while an 
additional source of lead, most probably, the s-process in AGB 



stars, is needed to explain the Pb abundances of the stars with 
[Fe/H] > -2.3. The non-LTE calculations for Th ii are presented 
in Sect.jSj Our conclusions are given in Sect. (6] 



2. Non-LTE calculations for Pbi 

To solve the coupled radiative transfer and statistical equilibrium 
(SE) equations, we use d a revi sed version of the DETAIL pro- 
gram (Butler & GiddingUm) based on the accelerated lambda 



iteration (AL I) method as described by 'Rybicki & H ummei 
[T992I) . The update was presented by IMashonkina et al 



1991 



2011 



). The non-LTE calculations were performed with the 



MARCS model atmospheres (iGustafsson et al.ll2008l PI. 



2.1. Pb I atomic structure calculations 

The atomic structure of neutral lead has been well studied with 
laboratory measurements. The ground level of Pb i is 6s~6p' ■'Pg, 
and the ionization of the 6p electron gives the ionic core 
(Pbn) 6s^6p, which splits into two levels, 'P^p (lower ion- 
ization limit 59819.57 cm"') and "P^^^-, (upper ionization limit 
73900.64 cm ') in the LS coupling scheme. The Pbi energy 
structure more closely approaches the J J than the LS coupling. 
Thus, the ground configuration 6s^6p^ possesses 5 energy lev- 
els which are designated as (1/2, l/2)o, (l/2,3/2)i, (3/2, 1/2)2, 
(3/2, 3/2)2, and (3/2, 3/2)o. The energy separation between the 
6s^6p^ levels is remarkably large, for example, the second and 
the fifth level have an excitation energy of /iexc = 0.97 eV and 
3.65 eV, respectively. Almost all the remaining known energy 
levels of Pbi belong to the 6s^6pnl electronic configurations. 
Hereafter, 6s- is omitted in the configuration designations. 

The early data on th e Pbi energy levels obtained before 
1958 were summarized by iMoord ( 1 1 95 8h . Afte r that, the analysis 
was g reatly extended u sing classical emiss ion dWood & Andrew! 
1 19681) and absorption dBrown et al.lll97'7h spectroscopy as well 
as selective laser excitation (the papers relevant to our study are 
cited below). The measurements of long Rydberg series termi- 
nating on the first ionizati on limit provided m any odd parity lev- 
els of the 6pns and 6pnd dBrown et al ] ll977l up to n ^ 59 and 
n - 7 7, respectively), an d 6png (n - 5-9 in IWood & Andrew! 
d 19681) and n = 21-47 in iDembczyriski et al.l dl994!) ) configu- 
rations. Even parity levels of the 6pnp (n = 7-60) and 6pnf (n 
- 5-56 ) configurations were measured by iHasegawa & Suziikil 
dT996l) a total of 205 energy levels with a total angular momen- 
tum J -0, 1, and 2. Information on several hundred autoionizing 
levels above the first ionization limit is available (see Ahad et al] 
I2OO5L and references therein). 

The d ata on transition probabilities in Pb i avai lable in the 
Hterature dBiemont et al.ll2000l iFuhr & Wies3ll992!) are incom- 
plete, and we rely on oscillator strengths computed in this study. 
The a tomic structu re of Pbi was calculated using the Cowan 
code dCowan!ll98l]) . The 22 even parity interacting configura- 
tions 6p^ + 6pnp (7 <n< 15) + 6pnf (5 <n< 12) + 6p6h -1- 6p7h 
+ 6s6p-7s + 6p'* were included in the energy matrix, while the 
odd set was restricted to 24 configurations, i.e., 6pns (7 <n< 
12) -I- 6pnd (6 <n< 12) + 6png (5 <n< 12) + 6sp-7p + 6s6p^ + 



5d^6s~6p^. The calcu lations generally follow the scheme used by 
iBiemont et al.l d2000l) . During the fitting of the calculated energy 
levels to the experimental ones, the number of varying energy 
parameters was restricted as much as possible. All the parame- 
ters describing the configurations in the 6pnl series for different 
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n were varied collectively keeping the ratios of the correspond- 
ing ab initio values. The exchange Slater integrals G* within the 
configurations were also varied in a similar way. The electro- 
static parameters including the configuration interaction param- 
eters, which were not optimized in the fitting procedure, were 
scaled down by a factor of 0.80, while the spin-orbit integrals 
were fixed at their ab initio values. As a result, the 67 even par- 
ity levels were described by the 25 parameters with a standard 
deviation of 73 cm ' and the 51 odd parity levels by 21 param- 
eters with a standard deviation of 74 cm"' . Our atomic structure 
calculations complemented the system of known energy levels 
in the region above 6.87 eV with the 20 even parity levels of the 
6pnf (n = 9-12, y = 3 and 4) and 6pnh (n - 6, 7) configurations 
and 29 odd parity levels belonging mostly to the 6png (n = 5-12) 
configurations. 

The wavefunctions obtained after the fitting of energy lev- 
els were used for the calculations of transition probabilities. 
The dipole transition integrals for all transitions were kept at 
ab initio Hartree-Fock values. The accuracy of our calcula- 
tions was estimated by the comparison with the experimen- 
tal data for common levels. Figure[T] shows that the differ- 
ence in l ifetime betwee n our calculations and the measure- 
ments of iBiemont et al.l (l2000l) is mostly within 75 %, with 
the only outlier being the measured to calculated lifetime ra- 
tio of about 2.3. On average, the predicted lifetimes approach 
to the measured values. It should be noted that the lifetimes 
of the highest 6pnd levels with J - 2 are critically dependent 
on a position of strongly interacting 6s6p^ level. It seems 
that the only level of the 6s6p'' configuration is firmly lo- 
cated at Eexc = 85870 crn'' JXrause et al.l [19861 iMuller et al.l 
II990I) . iDembczvhski et ai] (II994 estimated the position of the 
6s6p^ ^D" level at £exc = 68943 cm from analysis of the inter- 
actions in the 6pns and 6pnd series. They also established that 
"the state 6s6p^ does not occur as a dominant component in 
any of the electronic levels but it mixed into many odd levels". 
Its unperturbed value should be around 59300 cm To predict 
the levels of the 6s6p'' configuration, the same scaling of ab ini- 
tio values was applied to the parameters within the 6p'' subshell 
as that for the 6s^6p^ configuration. The average energy and the 
G'(6s,6p) parameter were adjusted to fit the energies of the ^Pj 



and D° levels. The energy of the 6s6p level appeared to 
be 58700 cm~', in quite g ood agreement with the estimates of 
iDembczvnski et all (Il994l) . Similar consistency between theory 
and measurements was found for the even parity 6pnp (7 <n< 
1 4) and 6pnf (6 < n< 8) levels when using the measured lifetimes 
of lLiet ail (120011) . 

Figure[T] shows also the me asured to calcula ted lifetime ra- 
tios from the calculations of Bie mont et al. I (l2000l) . who took into 
account the core polarization effect that is expected to give rise 
to the precision of the predicted lifetimes. As can be seen, our re- 
sults obtained without taking into accou nt the core polarizatio n 
agree very well with the calculations of iBiemont et al.l (l2000l) . 
Therefore, ^/-values calculated in this study are expected to 
have the same accuracy as the values in IBiemont et al.l ( l2000l) 
for all the transitions in Pb i. 

2.2. Lead model atom 

Energy levels. To construct the model atom, we used all the mea- 
sured and predicted levels in Pb i. They are shown in Fig.|2l The 
high-excitation (/Sexc ^ 6.04 eV) levels with common parity and 
energy separation of smaller than 0.01 eV were combined into a 
single level. The final model atom is fairly complete. It includes 




Fig. 1. Measured (IBiemont et al.ll2000l) to calculated Pbi life- 
time ratios for the 6pns (7 <n< 12) and 6pnd (6 <n< 13) levels, 
mostly with 7-1 and 2. Open ci rcles and crosses correspond to 
our and the IBiemont et al.l (l2000l) calculations, respectively. 
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Fig. 2. Measured and predicted energy levels of Pb i that form 
our final Pb model atom. 



97 levels of Pb i and the ground state bs^bp ^Pi^2 
second sublevel of the ground term, '^^^12 (^exc = 1 -75 e V), was 
taken into account only in the number conservation equation. 
The excited electronic configurations of Pb 11 produce levels with 
an excitation energy of more than 7.2 eV, and were thus ignored. 

Radiative rates. For 141 transitions in Pbi that connect the 
6p^ and 6p7p levels to the 6pns and 6pnd (n< 13) ones, we em- 
ployed accurate ^/-values based on measured n atural radiative 
lifetim es a nd theoretical branch ing ratios from IBiemont et al.l 
(I2OOOI) and lFuhr& wTes3(II992l) . For the remaining majority of 
transitions, we used oscillator strengths computed in this study. 
In total, 1 954 radiative bound-bound {b — b) transitions in Pb i 
were included in our SB calculations. 

For the Pb i ground state, we applied the photoioni zation 
cross-sections o-pj, measured bv lHeppinstall & Marrl(ll969l) from 
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the ionization threshold at 1671 to 1470 A. For the remaining 
levels, (Tph were computed using the hydrogen approximation 

formula with n - n^g, where ngff - Z -s/xn/Xi is the effective 
principal quantum number for the level with the ionization en- 
ergy Xi- Here, Z - 1 and xh is the hydrogen ionization energy. 
Such a choice was justified by the comparison of the experimen- 
tal threshold cross-section o-th^ e.^p = 10 Mb for the Pbi ground 
state with the cross-sections calculated with various n. The use 
of n = 6 results in the too low hydrogen approximation cross- 
section cTihj-ijyj = 0.0007 Mb. With neff = 1.35, we computed 
o"thr,/ivrf =1-3 Mb, which is much closer to the experimental data, 
although it is about a factor of 10 lower. 

Collisional rates. All levels in our model atom are coupled 
via collisional excitation and ionization by electrons and by neu- 
tral hydrogen atoms. The calculations of collisional rates rely on 
theoretical approxi mations. For electron-i mpact excitation, we 
used the formula of lvan RegemorteJ (Il962h for the allowed tran- 
sitions and we assumed that the effective collision strength Y - 
1 for the forbidden transitions. Electr on-impact ioniz ation cross- 
sections were computed according to lDrawi III (1196 Ih . 

For collisions with Hi atoms, w e employed the classi- 
cal D rawin formula as described by ISteenbock & Holwegeij 
(Il984l) . Over the past few decades, this formula has been crit- 
icized for significan tly overestimating the collision rates (see 
iBarklem et alj|201ll and references therein). At the same time, 
the need for a thermalizing process not involving electrons in 
the atmospheres of, in particular, very metal-poor stars, was in- 
dicated by many non-L TE spectral line formation studies (see 
iMashonkina et al.ll201 iL and references therein). Because no ac- 
curate calculations of either inelastic collisions of lead with neu- 
tral hydrogen atoms or other types of processes are available, 
we simulate an additional source of thermalization in the at- 
mospheres of cool stars by using parametrized Hi collisions. 
The Drawin formula can only be applied to allowed b - b and 
b - f transitions. For the forbidden transitions, a simple re- 
lation between hydrogen and electron collisio n rates , Ch - 
Ce ^J{mJmH)NHlNe, was used following Takedal (ll995l) . The SE 
calculations were performed with various efficiencies of colli- 
sions with H I atoms by applying a scaling factor of S'h = (no 
hydrogen collisions), 0.1, and 1. 

We discuss below the effects of the uncertainty in the adopted 
photoionization cross-sections and collisional rates on our final 
results. 

2.3. Departures from LTE for Pb i 

We first inspect the statistical equilibrium of lead in the two 
model atmospheres representing the solar atmosphere, with 
reff/logg/[Fe/H] = 5780/4.44/0, and the atmosphere of a typi- 
cal very metal-poor (VMP) cool giant (4500/1. 0/-3) from the 
calculations with Sh = 0.1. The departure coefficients, bj - 
^jNLTE^j^LTE^ for the Pbi and Pbii levels are shown in Fig. [3] 
Here, n^^^^ and n^^^ are the statistical equilibrium and TE num- 
ber densities, respectively. In these computations, the Pb abun- 
dance was adopted at logepb - 1.92 in the solar model and 
at logepb = -0.08 ([Pb/Fe] = -Hi) in the VMP model. In each 
model, lead is almost completely ionized, such that the fraction 
of Pb I nowhere exceeds 4% in the solar atmosphere and 0.6% 
in the cool VMP atmosphere. This explains why the non-LTE 
mechanisms for Pbi are similar in the two models. We summa- 
rize our findings as follows. 

1 . In the atmosphere outside log T5000 = +0.2 and +0.5 for the 
solar and VMP models, respectively, the most populated 6p~ lev- 
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Fig. 3. Departure coefficients, b, for the lowest 60 levels of Pb i 
and the ground state of Pb 11 as a function of log T5000 in the 
model atmospheres 5780/4.44/0 (top panel) and 4500/1. 0/-3 
(bottom panel). The first 21 levels are quoted in the right part of 
each panel. For all the higher excitation levels, their behaviour 
is very similar to that of the 188p to 217d levels, with b ^ I in 
the solar and VMP models inside log T5000 = -3 and -1, respec- 
tively. S'h = 0. 1 was used throughout. The vertical lines indicate 
the locations of line core formation depths for Pb 1 4057 (1) and 
3683 A(2). See text for more details. 



els of Pb I are underpopulated relative to their TE number densi- 
ties and they are close together. The main non-LTE mechanism 
is the ultraviolet (UV) overionization caused by superthermal ra- 
diation of a non-local origin below the thresholds at /l,/,^ - 1923, 
2033, and 2606 A for the three excited 6p^ levels. 

2. The populations of the 6p7s levels labeled in the model 
atom as 67s and 77s are controlled by the two mechanisms. 
The UV overionization (/},/„■ - 4022 and 4076 A) tends to drain 
the populations far inside the atmosphere. The competing pro- 
cess is radiative pumping of the strong 6p^ - 6p7s transitions 
(1/2, l/2)o - (1/2, l/2)i° (16p - 77s in our denotations, which 
produces the line at 2833 A), (l/2,3/2)i - (1/2, l/2)o° (26p - 
67s, 3683 A), and (l/2,3/2)i - (1/2, l/2)i° (26p - 77s, 3639 A) 
by the ultraviolet Jy-By{Te) excess radiation which tends to pro- 
duce enhanced excitation of the upper levels in the atmospheric 
layers, where the line wing optical depth drops below 1 . The net 
effect is that the 67s and 77s levels are underpopulated far inside 
the atmosphere, at log tsooo = -0.7 to +0.2 and -0.3 to +0.2 in 
the models 4500/1. 0/-3 and 5780/4.44/0, respectively, and they 
are overpopulated in the higher layers. No process seems to com- 
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pete with radiative pumping of the 26p - 157s (2476 A) and 36p 
- 157s (2663 A) transitions resulting in strong overpopulation of 
the upper level in the atmosphere outside log T50()() — +0.2. 

3. The overpopulation of the 6p7s levels is redistributed to 
the 6p7p levels through inelastic collisions. 

4. All the high-excitation levels with E^xc > 6.1 eV (starting 
from the 18th level in the model atom) couple thermally to the 
Pb n ground state. 

5. Pb n represents the state in which the majority of the ele- 
ment exists and it preserves, therefore, the TE population. 

The formation depth has been specified as the depth where 
TyAy" - 1, following iMihalasI (Il978h . Here, Ay is the ratio of 
the thermal to total absorption coefficient. Everywhere in the so- 
lar atmosphere. Ay ^ I, irrespective of the frequency point. In 
the cool VMP model. Ay < 1 due to significant contribution of 
Rayleigh scattering to the total opacity. For example, at the depth 
point log 75000 = 0, /ly - 0.7 at 4057 A and it decreases towards 
shorter wavelengths. Even lower Ay values are obtained in the 
higher atmospheric layers. 

The two lines, Pbi 4057 A (6p' (1/2,3/2)2 - 6p7s 
(1/2, 1 /2)i° denoted as 36p -77s in ourmodel atom) and 3683 A 
(6p2 (1/2, 3/2)i - 6p7s (1/2, l/2)o° denoted as 26p - 67s), and, 
in many cases, only one of them are used in stellar Pb abun- 
dance analyses. In each model, non-LTE leads to a weaken- 
ing of the Pb I lines compared to their LTE strengths and thus 
positive non-LTE abundance corrections. This is due to overi- 
onization (bi < 1) and also owing to the line source function 
5'/,/ - bu/bjBy rises {bu > bi) being above the Planck func- 
tion By in the line-formation layers (Fig. [3]). Here, b,, and bi are 
the departure coefficients of the upper and lower levels, respec- 
tively. We found that the departures from LTE are significantly 
larger in the VMP model compared to the solar one. The non- 
LTE abundance correction, Anlte = logfiNLTE - logSLXE, de- 
creases from +0.21 down to +0.09 dex for Pbi 3683 A in the 
solar model when Sh moves from to 1. Similar values, +0.19 
and +0. 1 1 dex, were obtained for Pb i 4057 A. In the 4500/1/-3 
model, ANLTE(Pb 1 4057 A) = +0.49 and +0.38 dex when 5h = 
and 1, respectively. 

Using 5h = 0.1, we computed Anlte for the grid of model 
atmospheres characteristic of the Galactic halo stars. Table [1] 
presents the resulting values for Pb i 4057 A. It can be seen that 
the non-LTE correction grows toward lower metallicity from 
Anlte - 0.16 dex in the solar model up to Anlte = 0.62 dex in 
the 5000/1. 50/-3 one. For a given metallicity, Anlte increases 
with increasing T^s and decreasing log g. 

Besides different prescriptions of the collisions with hydro- 
gen atoms, two different test calculations were performed for the 
4500/1/-3 model to assess the influence of crucial atomic data 
on final results. Here, we used Su- 0.1. First, we varied pho- 
toionization cross-sections by employing the principal quantum 
number n instead of Ueft in the hydrogen approximation formula. 
For the 6p~ and 6p7s levels, this resulted in the reduction of cr^,/, 
by a factor of 150 to 70. However, the effect on Anlte for Pbi 
4057 A was found to be minor, of less than 0.01 dex, because 
the departures from LTE for this line are controlled by the b - b 
radiative transitions between the 6p^ and 6p7s levels. 

Second, we checked the atomic model, where H i collisions 
were taken into account for the allowed transitions, but were 
neglected for the forbidden ones. As expected, this resulted in 
strengthening the departures from LTE. For Pb i 4057 A, Anlte 
grew from 0.41 dex up to 0.53 dex. Thus, poor knowledge of 
collisions with H i atoms is the main source of the uncertainty in 
the calculated non-LTE abundance corrections for the Pb i lines. 



Table 1. Non-LTE abundance corrections (dex) for the Pbi 
4057 A and Euii 4129 A lines from the calculations with 5h = 
0.1. 
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4750 
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-3 


0.53 


0.06 


5000 


1.5 


-3 


0.62 


0.07 



' Abundances used in non-LTE calculations. 



3. Solar lead abundance 

As a test and first application of the Pbi model atom, we de- 
termined the lead abundance for the Sun. We used solar central 
intensity observations taken from the Kitt Peak Solar Atlas of 
iBrault et al.1 (|1972|) . The only useful line of Pb i at 3683.46 A lies 
in a rather crowded spectral region, and the element abundance 
was derived applying the line-profile fitting method with the 
SIU program (lReetzll991 [). We employed the semi-empirical so- 
lar model atmosphere of [Holweger & Mueller ( 1974, hereafter, 
HM74) and the theoretical MARCS model dOustafsson et al.1 
,2008i) with Teif = 5780 K and \ogg = 4.44. A depth-independent 
microturbulence of 0.9 km s"' was adopted. 

The continuum level in the spectral range 3678-3688 A 
was fixed us i ng an atomic line list from the VALD database 
( Kupka et al.l 1 19991) and a molecular line list from iKurucd 
(Il994h . Our synthetic intensity profiles were convolved with a 
radial-tangential macroturbulence profile of Vmm: - 2.5 km s"'. 
Figure|4]illustrates the fit of the 3683 A blend, in particular, Pb i 
3683.46 A and the strong Fei 3683. 612 and 3683.629 A lines. 
For Pbi 3683.46 A, \oggf = -0.52 dBiemont et al.ll2000l) was 
employed throughout. 

Under the LTE assumption, we obtained logepbxTE = 1-96 
with the HM74 model and logepb,LTE = 1-85 with the solar 
MARCS model. It should be stressed that the solar photosphere 
lead LTE ab undance is lower than the meteoritic one log 
^ 2.0 4+0.03 dAsplund et al.l2009l) and 2.06+0.03 dLodders etal.l 
l2009h . independent of the applied model atmosphere. The two 
recent LTE determinations y ielded log gpb.LjE - 1 .75+0. 10 based 
on the theoretical 3D model f Asplund et ai.''2009') and 2.00+0.06 
based on the HM74 model (Lodders et al. 2009). 

We found that non-LTE decreases the discrepancy between 
the solar photosphere and meteoritic Pb abundance. When using 



5 



Mashonkina et al.: Non-LTE effects on stellar Pb and Th abundance determinations 



1.0 
0.9 

g 0.8 

B 

f> 0.7 



■S 0.6 
02 



0.5 
0.4 









j 






\ 

\ 

r \ 


A 

/ — 

1 : 

j. 

i 

i 

1 






Pbl ■ 

i 


i : 
/ 








■/ 






, , , I , , . 





3683.3 3683.4 3683.5 3683.6 
Wavelength (A) 



Fig. 4. Synthetic non-LTE {Sh - 0.1, dashed curve) and 
LTE (dash-dotted curve) disk-center intensity profiles of Pbi 
3683.46 A computed with the HM74 model atmosphere and 
loggpb = 2.09 co mpared with the observed spectrum of the 
iBrault et"aD (Il972h solar atlas (bold dots). Dotted curve corre- 
sponds to the atmosphere without lead. 

the HM74 model, the non-LTE abundances amount to log spb - 
2.12, 2.09, and 2.03 from the calculations with Sh - 0, 0.1, and 
1 , respectively. Figure|4]shows the best non-LTE fit for Sh -0.1. 
For the solar MARCS model, the non-LTE abundances are 2.00, 
1.97, and 1.92 depending on the used Sh value, i.e., for 0, 0.1, 
and 1, respectively. We note that, in all cases, the non-LTE cor- 
rection calculated for the emergent intensity spectrum is smaller 
than that for the flux spectrum. This is because the disk-center 
radiation emerges from deeper layers, where the departures from 
LTE are weaker From the analysis of the Pb i 3683 A line alone 
we can not decide about the efficiency of poorly known inelastic 
collisions with H i atoms in SE calculations, and we choose Sh - 
0. 1 based on our empirical esti mates for Ca i-Ca n and Fe i-Fe ii 
dMashonkina et al.ll2007l I20TTI) . Thus, our final solar lead abun- 
dance is logspb,NLTE - 2.09 with the HM74 model atmosphere 
and log epb.NLTE = 1 -97 with the solar MARCS model. 

4. Constraints on the Pb production mechanisms in 
the early Galaxy 

In this section, we revise the Pb abundances of metal-poor stars 
available in the literature and investigate the metallicity be- 
haviour of Pb/Eu abundance ratio. This ratio is particularly sen- 
sitive to whether the production of lead occurred in a pure r- 
process or whether there was a non-negligible s-process contri- 
bution. The r-nuclei of lead are produced through several decay 
channels. At the time when the r-process event stops (likely after 
a few seconds), there is the direct decay of very neutron-rich 
isobaric nuclei with A = 206 - 208 to ^o^Pb, 20Vpb, and ^o^Pb. 
Then there is the a- and yS-decay of nuclei with 210 < A < 231 
and A = 234 back to Pb, which operates shortly after the ter- 
mination of the r-process, within few million years. And fi- 
nally the radioactive decay of the long-lived Th and U isotopes 
produced Pb o ver the course o f almo st the age of the Galaxy. 
As stressed by iRoederer et al.l (l2009l) . no Pb isotopes can be 
measured individually due to the relatively small isotope shifts 
(< 0.02 - 0.03 A) and overall weakness of the Pb i lines. 

Hence, for the total elemental abun dances, the classi- 
cal waiting-point (WP) r-process model of iKratz et al.l (l2007h 
predicts loger_o(Pb/Eu) = 0.61 in a newly born star and 
log Sr, 13 (Pb/Eu) = 0.70 in a 13 Gyr old star (IRoederer et all 



I2009h . Another estimate of a pure r-process production ra- 
tio can be obtained using the Solar r-residuals. With an s- 
proce ss contribution to the solar P b of 81% (T ravag lio et al.l 
1200 llPI and to the solar Eu of 6 % dTravagho et al.ltl999r we 



computed the Solar System r-pro cess ratio loggssr(Pb/Eu) 
0.82. When computing s-fractions. lTravaglio et al.l (1 19991 



2001b 



used th e meteoritic isotope abundances from Anders & Grevessa 
(Il989l) . For comparison, the ratio of the meteoritic total abun- 
dances is logSmetCPb/Eu) = 1.51 and 1.53 acco rding to 
lAnders & Grevessel (119891) and iLodders et aP (l2009l) . respec- 
tively. 

The larg est data set on ste llar Pb abundances was com- 
piled by iRoederer et al.l (120101) . Their data are based on the 
analysis of Pbi 4057 A. Before correcting these abundances, 
we excluded all stars with only upper limits for Pb abundance, 
the carbon-enhanced stars with [C/Fe] > 0.3, and also VMP 
stars enriched in s-process material, with loge(Pb/Eu)> +1.8 
at [Fe/H] < - 2 (see s mall blue "x" in their Fig. 3) from the 
IRoederer et al.l (1201 Ol) sample. Regarding obtaining carbon 
abundan ces for these s tars to check upon, we also relied on the 
iSimmerer et al.l (l2004h study, which was used by Roede rer et al.l 
(120 lOl) to compile most of their stellar sample. In total, 49 
stars were selected in the -2.26 < [Fe/H] < -0.59 metallicity 
range . We added to t h em the two r-II stars, i.e., CS 31082- 
001 ('Hill et al.' '2 0021 iPlez et al.l |2004 and HE 1523-0901 
(iFrebel et al.i.200 7 ), 2012, in preparation) , and a lso the halo star 
HD 29907 from lSitnova & Mashonkiiial (1201 ll) . The stars used 
are listed in Table |2| with stellar parameters and LTE element 
abundances taken from the cited papers and with non-LTE 
abundance corrections obtained in this study. The non-LTE cal- 
culations were performe d with the MARCS model atmospheres 
(iGustafsson et al.ll2008l) . Where necessary, the MARCS model 
structures were interpolated for given Teff, logg, and [Fe/H] 
using a FORTRAN-based routine written by Thomas Masseron 
(http : //marcs . astr o . uu . se/ sof tware . php). For each 
star, the SE calculations were performed using its given LTE 
element abundance. 

For self-consistency, we also revised the europium abun- 
dances of the stars in our sample. The non-LTE calcula- 
tions for Euii were p e rforme d using the method treated by 
iMashonkina & GehrenI (l2000l) . Everywhere, Sh - 0.1 was 
adopted. In the stellar parameter range with which we are con- 
cerned, non-LTE leads to weakened Eu 11 lines and positive abun- 
dance corrections. We computed Anlte for the seven lines, i.e., 
Eu II 3724, 3819, 4129, and 4205 A with £exc = and Eu 11 3907, 
4435, and 4522 A with Eexc = 0.21 eV. Table [Upresents the re- 
sulting values for Eu 11 4129 A in the grid of model atmospheres 
characteristic of the Galactic halo stars. Table |2| does the same 
for the stars in our sample. For each line, Anlte grows toward 
lower metallicity and lower surface gravity, although, it is over- 
all small and does not exceed 0.17 dex for different lines. For 
example, Anlte ranges between 0.02 and 0.06 dex for different 
lines in the 5600/3.75/-0.59 model, between 0.04 and 0.14 dex 
in the 4825/1. 5/-2.9 model, and between 0.05 and 0.17 dex in 
the 4050/0. 0/-2. 08 model. For each model atmosphere, Anlte 
of Euii 4129 A turned out to well represent the mean non-LTE 
correction from all the computed lines. For a given star, we used 
it as the mean non-LTE correction and applied it to the observed 
europium LTE abundance. 



^ These authors actually indicated 91 %. However, using s-fractions 
of 90, 60, 77, and 89 % for the ^o-tpb, ^o^Pb, ^o'p^,^ 208p,, isotopes, 
respectively, from their Table 3, we calculated an s-fraction of 81 % for 
the total lead abundance. 
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Fig. 5. Stellar non-LTE (filled circles) and LTE (open circles 
for the stars with \ogg < 3 and triangles for log^ > 3) Pb/Eu 
abundance ratios as a function of metallicity. The solid line in- 
dicates the Solar system meteoritic value. The two dotted lines 
correspond to a pure r-process production of Pb and Eu predicted 
by the waiting-point r-process model for a newly born (the low- 
lyi ng curve) and 1 3 Gyr old (the up-lying curve) star as given 
bv lRoederer et alj (|2009|) . The non-LTE abundances were calcu- 
lated with 5h = 0.1. 



The non-LTE and LTE Pb/Eu abun dance ratios are plotted 
in Fig.|5l As shown bv lHill et al.l dlool for CS 31082-001 and 
iFrebel et al.l ( l2007h for HE 1523-0901, the heavy elements be- 
yond the iron group in these stars originate from the r-process, 
and the observed element abundance ratios can be used to test 
theoretical nucleosynthesis models. We found that both stars 
have very similar Pb/Eu abundance ratios at log eNLTE(Pb/Eu) 
= 0.67 and 0.68, respectively, and that th e WP r-process mode l 
for a 13 Gyr old star ( iKratz etal.i,2007i iRoederer et alj|2009l) 
perfectly reproduces the observations. We also stress that the 
Pb/Eu abundance ratios of the r-II stars match the Solar r-process 
rat io log £ssr ( Pb/Eu ) = 0.82 within the error bars. According 
to iPlez et ai] (l2004l) . the uncertainty in the Pb abundance of 
CS 31082-001 is 0.15 dex. This provides strong evidence for 
universal Eu and Pb relative r-process yields during the Galactic 
history. 

Despite the large scatter in the abundances, a clear upward 
trend is seen in the loge(Pb/Eu) - [Fe/H] plane for all stars 
with [Fe/H] > -2.3. The Pb/Eu abundance ratio grows from 
logeNLTE(Pb/Eu) ^ 0.8 at [F e/H] - -2 up to the s olar value 
at [Fe/H] = -0.6. Following iRoederer et all (1201 Ol) . we sep- 
arated these stars in two groups depending on their metallic- 
ity. The low-metallicity sample (hereafter, LMS) includes 29 
stars in the -2.26 < [Fe/H] < -1.40 metallicity range, and the 
mildly metal-poor sample (-1.4 < [Fe/H] < -0.59, hereafter, 
MMPS) includes 21 stars. It is worth noting that each subsample 
is highly surface gravity biased due to a common selection ef- 
fect. Namely, most LMS stars are cool giants, while the MMPS 
stars are mainly dwarf stars. 

As seen in Fig.|5j the LMS stars have significantly higher 
Pb/Eu ratios compared to that of the r-II stars, independent 
of either LTE or non-LTE. Excluding a clear outlier, i.e., the 
star HD 108317 ([Fe/H] = -2.18, logeNLTE(Pb/Eu) = 1.91), 
we calculated the mean logeLTE(Pb/Eu) - 0.93 ± 0.23 and 
logeNLTE(Pb/Eu) = 1.19 ±0.24 from the 28 LMS stars. We shall 
conclude that the s-process in AGB stars started to produce lead 
before the Galactic metallicity grew to [Fe/H] = -2.3. The low- 



metal l icity s-process nucle synthesis calculations (IRoederer et alJ 
120 Id . 'Lugaro et al .l 1201 1|) do not contradict this observational 
finding. Roederer et alj ( l2()Tol) presented the surface composi- 
tion at the end of the AGB phase for a number of intermediate- 
mass models with [Fe/H] = -2.3 and M = 4.5 - 6Mq. The min- 
imum value produced by those models is loge(Pb/Eu) - 1.9. 
Given sufficient time, their contributions could result in raising 
the Galactic Pb/Eu abundance ratio above the pure r-process one. 
Earlier, from analyses of t he ste llar ( Ba, La, Nd)/(Eu, Dy) a bun- 
dance ratios. iBurris et al.l (l2000l) and lSimmerer et al.l (|2004 ob- 
tained evidences for that the s-process might be active as early 
as [Fe/H] - -2.3 and -2.6, respectively. 

For the MMPS stars, the obtained mean, log eNLTE(Pb/Eu) = 
1.58 ± 0.31, is very close to the Sol ar System value, 
loge,net(Pb/Eu) = 1.53 ('Lodders et al. 2009^. This supports the 
theoretical result of iTravaglio et al. (200D who predicted that 
the AGB stars with [Fe/H] ^ -1 made the greatest contribu- 
tion to the solar abundance of s-nuclei of lead. It is worth not- 
ing that LTE leads to 0.22 dex lower Pb/Eu abundance ratio for 
the MMPS stars. The star G 58-25, with [Fe/H] = -1.40 and 
log eNLTE(Pb/Eu) = 2.28, was not included in the mean calcula- 
tions. 

Generally, our findings are stable with respect to a variation 
in the stellar parameters that coiTespond to reasonable observa- 
tional uncertainties and a variation in the Sh value. For example, 
fo r CS 3 1082-001, a change of -i-lOO K in T^ff (a ~ 2o- error 
in lHill et al. (2002)) produces a change in the calculated Pb non- 
LTE abundance of -i-0. 16 dex. Pb 1 4057 A is insensitive to a vari- 
ation in log g and micro turbulence velocity, Vmic- The total effect 
of varying T^ff, log g (by -1-0.3 dex), and Vmic (by -H0.2kms~') on 
the derived Eu abundance was estimated bv lHill et alj (|2002|) as 
to be -1-0. 12 dex. Thus, the Pb/Eu abundance ratio is only slightly 
affected by the uncertainties in stellar parameters. The non-LTE 
calculations performed for CS 31082-001 with = and 1 
result in only 0.04 and 0.02 dex lower log e(Pb/Eu) values com- 
pared to that for Sh - 0.1. To evaluate the uncertainty in the 
Pb/Eu abundance ratios of the LMS and MMPS stars, we cho- 
sen the two representative star s HD 3008 (4250/0.25/ -2.08) and 
G 102-020 (5250/4.44/- 1 .25).|SiminerereLalJ dlool estimated 
that the stellar parameters of their sample are internally consis- 
tent to AFeff = 100 K, Alog g = 0.25, and AVmic = 0.1 kms^'. 
An upward revision of the stellar parameters for HD 3008 results 
in a cumulative change of +0. 1 1 dex in log epb and H-0. 10 dex in 
log eeu- For G 102-020, the abundance errors caused by the un- 
certainties in the stellar parameters also have a common sign for 
Pb and Eu and amount to 0.09 and 0.06 dex, respectively. 

We note that we have not considered effects on the abun- 
dances as a result of using 3D hydrodynamical model atmo- 
spheres. However, those effects were estimated for Pbi 4057 A 
(Eexc - 1.32 eV) and the Eun Unes arising from the ground 
sta te using the 3D- ID LT E abundance corre ctions calculated 
bv lCoUet etaP (120071) and iHavek etaP (l201ll) for a sample of 
"fictitious" atomic (Na i, Mg i, Ca i, Fe i, and Fe ii) lines at se- 
lected wavelengths, with varying lower-level excitation poten- 
tials and line strengths. Overall, 3D-1D corrections are negative 
for the lines of neutral atoms and decrease in absolute value with 
increasing fiexc- In the 4717/2. 2/-1 models, they have similar 
values for all the chemical species, i.e., 3D-1D varies between 
-0.20 and -0.28 dex for the lines with £exc = and between 
-0.10 and -0.18 dex for ^exc = 2eV. The 3D effects grow to- 
ward lower metallicity, however, with different rate for differ- 
ent atoms. In the 4732/2.2/-2 and 4858/2.2/-3 models, the 3D- 
ID corrections are similar for Mgi and Fei, but different from 
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that for Nai and Cai. For the weak Nai line with iSexc = 1.3 eV 
and an equivalent width of EW < 20 mA, 3D- ID ^ -0.12 dex 
in the -3 < [Fe/H] < -1 metallicity range. For the weak Fei 
line with the same excitation potential, 3D-1D ^ -0.4, -0.3, and 
-0.17 dex in the [Fe/H] = -3, -2, and -1 model, respectively. 
The 3D effects are overall weak for the Fe n lines arising from the 
ground state. For example, for the moderate strength line (EW 
^ 50 mA), 3D-1D = -0.03, +0.15, and +0.15 dex in the [Fe/H] 
= -3,-2, and -1 model, respectively. We assume that 3D-1D 
corrections for the investigated Eu n lines are similar to that for 
the "fictitious" Fe ii line with /iexc = eV. The 3D- ID correction 
for Pb I 4057 A lies, probably, between those for Na i and Fe i. 
In case it is as large as that for Na i, accounting for the 3D ef- 
fects results in approximately 0.1 dex lower Pb/Eu abundance 
ratio for the two r-II stars and 0.25 dex lower Pb/Eu value for 
the low-metallicity and mildly metal-poor samples. With such 
changes, our conclusions remain valid. If the 3D- ID correction 
for Pb 1 4057 A is as large as that for Fe i, the mean Pb/Eu abun- 
dance ratios of the r-II stars and the LMS stars are expected to be 
approximately 0.4 dex lower compared with the corresponding 
NLTE+ID values. However, in this case, it would be hard to un- 
derstand why the Pb/Eu abundance ratio of the strongly r-process 
enhanced stars is significantly lower (by 0.5 dex) compared with 
the Solar System r-process value. 

We note here that full 3D non-LTE computations have so 
far only been performed for Li i in the Sun a nd metal-poor stars 
dAsplund et al'l l2003l iBar klem et all I2003L ICavrel et al.1 l2007l 



Sbordon e et al.l 20101) and for Oi in the Sun ( Asplund et al.l 
2004). It was found that the 3D+non-LTE effects are overall 
small for the Li i resonance line because the line strengthening 
due to the cooler temperatures in the upper atmospheric layers 
of the 3D model on the one hand, and line weakening from in- 
creased overionization on the other hand, largely cancel each 
other out. For O i, the non-LTE abundance corrections turned out 
similar in the 3D and I D solar model atmospheres. In contrast, 
IShchukina et al.l (l2005h predicted stronger non-LTE effects for 
Fei in 3D than in the ID model, based on their 1.5D+non-LTE 
calculations for the Sun and HD 140283 (5700/3.7/-2.5). Full 
3D non-LTE computations need to be done for Pb i to decide 
what is a difference in Pb/Eu abundance ratio between 3D+non- 
LTE and ID+non-LTE cases and how this might affect our con- 
clusions. 

5. Non-LTE calculations for Th ii 

5.1. Thorium model atom 

Energy levels. In the range of stellar parameters that we are con- 
sidering here, the majority of the element exists as Th ii. For ex- 
ample, the fraction of neutral thorium, with its ionization energy 
Xion = 6.3 eV, does not exceed lO"-' throughout the 4500/1/-3 at- 
mosphere and 610"^ throughout the solar atmosphere. We thus 
neglected Th i when constructing the model atom. The Th ii term 
structure is produced by multiple electronic configurations (see 
Fig. |6ll and consists of thousands energy levels. The ground 
configuration is 6d~7s. We used the energy levels with /iexg 
< 6.88 eV, which amount to 416 levels from Blaise & Wyar| 
(Il992h . The higher excitation levels play a minor role in the pop- 
ulation and depopulation of Thii, because the next ionization 
stage, Thiii, represents a minor fraction of the thorium abun- 
dance. For example, the contribution of Thiii to thorium abun- 
dance is smaller than 0.2 / 1.5 % in the cool VMP / solar mod- 
els outside log rsooo = 0. The levels with common parity and 
close energies were combined whenever the energy separation 



is smaller than 0.01 eV at /Sgxc < 4 eV and smaller than 0.1 eV 
at higher ZSexc. The final model atom includes 184 levels of Thii 
(Fig.©. 

The Thm ground state 5f6d^H4 is a poor representation 
of the Thm ionization stage. The odd 5f6d and 5f7s and the 
even 6d- and 6d7s configurations of Thm produce many low- 
excitation levels. The energy of 6(f ^¥2, fiexc - 63.2 cm"', is 
the smallest distance between the lowest levels of both parities 
ever found in an atomic spectrum. As a result, the Th m partition 
function (U) is significantly larger compared with the statisti- 
cal weight g - 9 of the ground state. For example, t/(Th m) = 
49.14 at the temperature T - 4990 K and grows to f/(Thm) = 
81.11 at T - 8318 K. These values we r e com puted using the 
Th m energy levels from lBlaise & Wvarll (Il992h . Assuming that 
the low-excitation levels of Thm are thermally coupled to the 
ground state, we included Thm in the model atom as a single 
state with g - 49, while keeping the energy of the ground state. 

Radiative rates. Accurate g/- values based on measured nat- 
ural radiative lifetimes and branching ratios were obtained by 
iNilsson et al.l (l2002h for 180 transitions between the lowest (^exc 
= - 2.25 eV) even parity levels of the 6d^7s, 6d7s^, and 6d^ 
electronic configurations and the intermediate energy (fiexc = 
2.26 - 4.08 eV) odd p arity levels. For anothe r 1244 transitions, 
we used g/- values of iKurucz & Belli (11995) whic h are acces- 
sible via the VALD database jKupka et al.l Il999l) . These two 
sources provide the data for 955 b - b transitions between com- 
bined levels in our model atom. They all were included in SE 
calculations. We refer to this atomic model as Th-1. 

The second atomic model, Th-2, includes the extended sys- 
tem of b - b radiative transitions. We estimated the total num- 
ber of allowed transitions in the system of 416 Thii levels as- 
suming that every transition between the different parity lev- 
els with Ay = 0, +1 is an allowed one. The exception are the 
J - ^ J - transitions. This gives 20214 transitions, and 
for 18 790 of them, oscillator strengths (fij) are missing. The 
1 562 of the latter transitions, with the energy separation of less 
than 0.25 eV, were treated as forbidden ones because of insignif- 
icant contribution of the radiative to the total transition rate, ir- 
respective of the transition ^/-value. For every of the remain- 
ing 17 228 transitions, fij = lO"'* was adopted. Such a choice 
was based on the insp ection of the known ;?/- values for Th n. 
No tra nsition from the iNilsson et"aD (l2002l) and iKurucz & Belli 
(Il995h data sets has fij < 10 In this atomic model, the 5 661 
b - b transitions were included in SE calculations. 

The photoionization cross-sections were computed using the 
hydrogen approximation formula with n = neff , because no accu- 
rate data is available for the Th 11 levels. 

Collisional rates. The same formulas as for Pb i were em- 
ployed for the transitions in Th n, with T = 2 for the forbidden 



5.2. Departures from LTE for Th 11 

The non-LTE calculations for Thn were performed with the 
MARCS solar model atmosphere and with the models repre- 
senting the atmospheres of VMP stars with [Fe/H] = -2 and 
-3 (see TableO. Figure|7]shows the departure coefficients of the 
selected Thii levels in the 5780/4.44/0 and 4500/1. 0/-3 models. 
We found that the lowest 72 levels, i.e., up to the level denoted 
as 726d (iSexc = 2.46 eV), keep nearly TE populations (b ^ 1) 
in the atmosphere inside log rsooo = -2 and - 1 for the solar 
and VMP models, respectively. For all the higher excitation lev- 
els, b > I. This is because the population of the odd parity 
5f6d^ and 6d7s7p levels with ^exc = 3 - 4.8 eV is controlled 
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> 
LU 



■ 1825f 



■ 1655f 



1495f* 



1345f 



: 1196d* 
- 986d* 



1235f* 



965f* = 



■ 1075f 



876d 



: 805f* 



= 726d 



-326d 



^ 505f* 



= 195f* 



M16d 



= 26d 



Th II 



-55r 



6d^ 6d' 6d' 6d7s 6d7s 5f6d 5f6d 5f6d 5f 5f 5f7s 5f7s 7s' 
7s 7p 7p 7s 7s 7p 6d 7s 6d 7s 7p 7p 



Fig. 6. Model atom of Th ii. The names of the levels important for understanding the non-LTE mechanisms are indicated. 



by the strong radiative transitions from the three lowest 6d-7s 
states and the 6d7s^ level, such as 6d-7s'*F3/2 - 6d7s7p (/Sexc - 
3.084 eV, 7 = 5/2) denoted in our model atom as 16d - 986d* 
and producing the line at 4019 A, 16d - 965f* (4094 A), 26d - 
1196d* (3741 A), and 36d - 1495f* (2895 A). These transitions 
are pumped by the ultraviolet Jy - By{Tg) excess radiation and 
produce enhanced excitation of the upper levels far inside the at- 
mosphere. This overpopulation is redistributed to the even parity 
high excitation levels through inelastic collisions. 

The Th-2 atomic model includes more pumping transitions 
compared to the Th-1 one, and the departures from LTE are ex- 
pected to be stronger when using the Th-2 model. We checked 
the populations of the lower and upper levels of the Th ii 4019 A 
transition. In the line formation layers, around log T5000 = -0.45 
and -0.55 in the solar and 4500/1. 0/-3 models, respectively, the 
difference in non-LTE populations between the Th-2 and Th-1 
atomic models is negligible for the Thii ground state. This is 
easy to understand because the lowest Thii levels contain the 
majority of the element, and no mechanism is able to signifi- 
cantly change their populations compared to the TE ones. For 
the upper level, 986d*, a small difference of 3 % between using 
Th-2 and Th-1 was found in the VMP model, and very similar 
populations were obtained in the solar model. 

The Thii lines used in abundance analyses all arise from 
either the ground state or low-excitation levels with ^exc ^ 



0. 51 eV, for which b ^ I holds in the line formation layers 
(Fig-0- Nevertheless, for each line, non-LTE leads to its weak- 
ening compared to the LTE strength, owing to the overpopu- 
lation of the upper level relative to the TE population that re- 
sults in b,i/bi > 1 and the rise of the line source function above 
the Planck function in the line-formation layers. Overall, the ob- 
tained non-LTE abundance corrections are positive. For the Sun, 
the departures from LTE are small irrespective of the applied 
Sn value, such that ANLTE(Thn 4019 A) = 0.06 and 0.00 dex in 
case of S'h = and 1, respectively. In the VMP models, the non- 
LTE correction is very sensitive to treatment of inelastic colli- 
sions with H I atoms. For example, for 4500/1/-3, it amounts to 
Anlte - +0.52 dex in the calculations with pure electronic colli- 
sions and it varies between Anlte = +0.12 and +0.03 dex when 
collisions with H i atoms are taken into account with Sh =0.1 to 

1. The use of the Th-1 atomic model results in 0.01 dex smaller 
non-LTE corrections. 

Table |3] presents Anlte for the two lines, Thii 4019 and 
4086 A, from the calculations with the Th-2 atomic model and 
Sh = 0.1. For a given Teff , the departures from LTE grow toward 
lower surface gravity due to decreasing contribution of the colh- 
sional to the total transition rates. Maximum Anlte were calcu- 
lated for the 4000/0.0/-2 model, and they amount to +0.21 dex 
for Th II 4019 A and +0. 14 dex for Th 11 4086 A. The metalUcity 
effect is less pronounced in Table[3]because, for the [Fe/H] = -3 
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Table 3. Non-LTE abundance corrections (dex) for the Thii 
lines from the calculations with 5h = 0. 1 . 



Fig. 7. The same as in Fig. [3] for the selected levels of Thii and 
the Th m ground state from the calculations with the Th-2 atomic 
model. The vertical lines indicate the locations of line core for- 
mation depths for Thii 4019. 



models, the non-LTE calculations were performed with a higher 
Th abundance compared with that in the [Fe/H] = -2 models. 

6. Conclusions 

We built a comprehensive model atom for neutral lead using 
atomic data for the energy levels and transition probabilities 
from laboratory measurements and theoretical predictions. The 
non-LTE calculations for Pb i were performed for the first time 
for the Sun and for a set of stellar parameters characteristic of 
metal-poor stars in the -2.95 < [Fe/H] < -0.59 metallicity 
range. We found that the main non-LTE mechanism for Pb i is 
the ultraviolet overionization, although radiative pumping of the 
strong transitions arising from the low-lying levels tends to pro- 
duce enhanced excitation of the Pbi levels with E^xc > 4 eV. 
Overall, non-LTE leads to weakened Pb i lines and positive non- 
LTE abundance corrections. The departures from LTE grow with 
decreasing metallicity and, for a given metallicity, increase to- 
ward higher effective temperatures and lower surface gravities. 
The main source of the uncertainty in the calculated non-LTE 
abundance corrections are poorly known inelastic collisions with 
Hi atoms. From analysis of either solar or stellar Pbi lines we 
can not empirically constrain their efficiency in SE calculations. 
With Sh- 0.1, chosen based on our estimates for Cai-Caii and 
Fe i-Fe ii ( Mashonki na et al. 2007, ,2011.) . Anlte of Pb i 4057 A 
ranges between 0.16 dex in the solar model and 0.56 dex in the 
4825/1. 5/-2.9 model. 
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' Thorium abundance used in non-LTE calculations 

^ Absolute thorium abundance, log eth 

' Calculated equivalent width is smaller than 1 mA 



We also built a model atom for singly-ionized thorium using 
atomic data for the energy levels from laboratory measurements. 
The non-LTE calculations for Th ii were performed for the first 
time for the Sun and for the small grid of model atmospheres 
with [Fe/H] — -2 and -3. In contrast to Pb i, Thii is the major- 
ity species in the stellar parameter range that we have covered, 
and the main non-LTE mechanism for Th ii is connected to the 
pumping transitions arising from the low-excitation levels, with 
Eexc < 1 eV. Overall, non-LTE leads to weakened Th n lines and 
positive non-LTE abundance coiTections. Compared with Pbi, 
the departures from LTE for Th ii are even more sensitive to the 
efficiency of collisions with Hi atoms. With S'h = 0.1, Anlte of 
the investigated Thn Unes nowhere exceeds 0.21 dex. 

We caution against using the LTE assumption in stellar anal- 
yses of Pb abundance and also the elemental abundance ratios 
involving Pb and some other element observed in the lines of 
their majority species, such as Eu, Th, or U. 

Taking advantage of our SE approach for Pbi, we ob- 
tained good agreement, within 0.03 dex, between meteoritic 
and absolute solar abundan ces from the calculations with the 
iHolweger & Muelled (11974 solar model atmosphere and S\i - 
0.1. Our final value is logepb,© - 2.09. 

We the n calculated th e non- LTE abundances of Pb and 
Eu for the iRoederer et alJ (1201 Oh metal-po or stars and com- 
plemented that s ample with HE 1523-0901 dFrebel et al.ll2007[) 
and HD 29907 dSitnova & Mashonkinall201 ih . We found that 
the two strongly r-process enhanced (r-II) stars have very sim- 
ilar Pb/Eu abundance ratios, with the mean log eNLTE(Pb/Eu) = 
0. 68+0.01, and the wait ing-point r-process model as presented 
bv lRoederer et al.l (l2009l) reproduces the observations very well. 
The revised Pb/Eu abundance ratios of the r-11 stars match, 
within the error bars, the corresponding solar r-process ratio. 
Thus, the universality of the r-process was proved not only for 
the secon d r-process p e ak elements from Ba to Hf as found 
earlier by ISneden et al.l (Il996l) . iHill et al.1 (l2002l) . ISneden et al.l 
(I2008L and references therein), but also for the heavier element 
Pb. 

It was found that the stars in the -2.3 < [Fe/H] < -0.6 
metallicity range, reveal a clear upward tre n d in th e log e(Pb /Eu) 
- [Fe/H] plane. Following iRoederer et al.l (1201 Ol) . we separated 
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these stars in two groups depending on their metallicity. The 
lower metallicity stars ([Fe/H] < -1.4) turned out to have, on 
average, 0.51 dex higher Pb/Eu ratios compared with that of the 
r-II stars. This led us to conclude that the s-process production 
of lead started before the Galactic metallici ty ^rew to [F e /H] — 
-2. 3. It is worth noting th at the studies of iBurris et al.l (|2000|) 
and lSimmerer et al.l (|2004 who obtained (Ba, La, Nd)/(Eu, Dy) 
abundance ratios of metal-poor stars, also found that the s- 
process has been active as early as [Fe/H] - -2.3 and -2.6, 
respectively. 

The mean, log eNLTE(Pb/Eu) - 1.58 + 0.31, for the mildly 
metal-poor subsample (-1.4 < [Fe/H] < -0.59) turned out 
very close to the So lar System value, logemet(Pb/Eu) — 1.53 
dLodders et al.ll2009h. T his is in agreement with the theoretical 
result of Travaglio et al.l (1200 Ih who predicted that the AGE stars 
with [Fe/H] ^ -1 made the greatest contribution to the abun- 
dance of s-nuclei of lead in the presolar nebula from which the 
Sun originated. 
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Table 2. Non-LT E and LTE abundance s of Pb and Eu and non-LTE abundance corrections (dex) for the Pb i 4057 A and Eu ii 
4129 A lines in the lRoederer et al.l (1201 Ol) stellar sample. 
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0.06 


HD 29574 


4250 


0.80 


-2.00 


0.25 


-0.06 


0.31 


-0.57 


-0.63 


0.06 


HD 29907^ 


5500 


4.64 


-1.55 


1.03 


0.70 


0.33 


-0.45 


-0.50 


0.05 


HD 37828 


4350 


1.50 


-1.62 


1.14 


0.77 


0.37 


-0.48 


-0.53 


0.05 


HD 44007 


4850 


2.00 


-1.72 


0.67 


0.31 


0.36 


-0.87 


-0.94 


0.07 


HD 63791 


4675 


2.00 


-1.90 


0.58 


0.22 


0.36 


-0.85 


-0.92 


0.07 


HD 74462 


4700 


2.00 


-1.52 


0.79 


0.49 


0.30 


-0.32 


-0.39 


0.07 


HD 105755 


5700 


3.82 


-0.83 


1.71 


1.43 


0.28 


0.07 


0.02 


0.05 


HD 106516 


6170 


4.21 


-0.81 


1.85 


1.56 


0.29 


0.05 


-0.04 


0.09 


HD 108317 


5234 


2.68 


-2.18 


0.69 


0.17 


0.52 


-1.22 


-1.32 


0.10 


HD 121135 


4934 


1.91 


-1.54 


0.72 


0.38 


0.34 


-0.63 


-0.70 


0.07 


HD 122956 


4508 


1.55 


-1.95 


0.19 


-0.13 


0.32 


-0.73 


-0.79 


0.06 


HD 141531 


4360 


1.14 


-1.79 


0.52 


0.20 


0.32 


-0.66 


-0.72 


0.06 


HD 166161 


5350 


2.56 


-1.23 


1.19 


0.84 


0.35 


-0.41 


-0.48 


0.07 


HD 171496 


4950 


2.37 


-0.67 


1.60 


141 


0.19 


0.17 


0.11 


0.06 


HD 175305 


4770 


1.80 


-1.73 


0.06 


-0.28 


0.34 


-0.82 


-0.89 


0.07 


HD 187111 


4270 


1.05 


-1.97 


0.21 


-0.10 


0.31 


-0.82 


-0.88 


0.06 


HD 201891 


5910 


4.19 


-1.09 


1.58 


1.25 


0.33 


-0.17 


-0.22 


0.05 


HD 204543 


4670 


1.49 


-1.87 


0.41 


0.05 


0.36 


-0.98 


-1.05 


0.07 


HD 206739 


4650 


1.78 


-1.72 


0.69 


0.38 


0.31 


-0.55 


-0.62 


0.07 


HD 210295 


4750 


2.50 


-1.46 


1.04 


0.72 


0.32 


-0.27 


-0.34 


0.07 


HD 214925 


4050 


0.00 


-2.08 


-0.13 


-0.50 


0.37 


-0.97 


-1.09 


0.12 


HD 220838 


4300 


0.60 


-1.80 


0.39 


0.05 


0.34 


-0.86 


-0.93 


0.07 


HD 221170 


4510 


1.00 


-2.16 


0.21 


-0.09 


0.30 


-0.80 


-0.86 


0.06 


HD 235766 


4650 


1.20 


-1.93 


0.42 


0.10 


0.32 


-0.80 


-0.86 


0.06 


HE 1523-0901 3 


4630 


1.00 


-2.95 


0.11 


-0.34-* 


0.45 


-0.57 


-0.62 


0.05 



' from the calculations with 5h = 0. 1 , 

^ stellar parameters and LTE abundances from Sitnova & Mashonkina (2011 
^ stellar parameters and europium LTE abundance from Frebel et al. (2007), 
* Frebel et al. (2012, in preparation). 
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